Topological properties of a time-integrated activity driven network 
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Here we consider the topological properties of the integrated networks emerging from the activity 
driven model [Perra at al. Sci. Rep. 2, 469 (2012)], a temporal network model recently proposed 
to explain the power-law degree distribution empirically observed in many real social networks. By 
means of a mapping to a hidden variables network model, we provide analytical expressions for the 
main topological properties of the integrated network, depending on the integration time and the 
distribution of activity potential characterizing the model. The expressions obtained, exacts in some 
cases, the results of controlled asymptotic expansions in others, are confirmed by means of extensive 
numerical simulations. Our analytical approach, which highlights the differences of the model with 
respect to the empirical observations made in real social networks, can be easily extended to deal 
with improved, more realistic modifications of the activity driven network paradigm. 

PACS numbers: 05.40.Fb, 89.75.Hc, 89.75.-k 



I. INTRODUCTION 

Modern network science allows to represent and ratio- 
nalize the properties and behavior of complex systems 
that can be represented in terms of a graph [H [5] . Re- 
search in this area has focused in a twofold objective: A 
data-driven effort to characterize the topological prop- 
erties of real networks [SHS], and a posterior modeling 
effort, aimed at understanding the microscopic mecha- 
nisms yielding the observed topological properties PQ [3] , 
as well as the effects that a complex topology has on 
dynamical processes running on top of it [H [7] . 

Until recently, a large majority of work in the field 
of network science has been concerned with the study 
of static networks, i.e. networks in which topological 
properties do not change in time. Presently, however, a 
lot of attention is being devoted to the temporal dimen- 
sion of networked systems [5]. Indeed, many real net- 
works are actually dynamical structures, in which edges 
appear, vanish or are rewired at different times scales. 
An important example is given by social networks [9] , in 
which social relationships are represented by a succession 
of contact or communication events, continuously created 
or terminated between pairs of individuals. In this sense, 
the social networks previously considered in the literature 
[10H12J represent an projection or temporal integration of 
time- varying graphs, in which all the links that have ap- 
peared at least once in a time integration window T are 
present in the projection. 

The recent availability of large digital databases and 
the deployment of new experimental infrastructures have 
made possible the real-time tracking of social interac- 
tions in groups of individuals and the reconstruction of 
the corresponding temporal networks |13H16j . The newly 
gathered empirical data poses new fundamental questions 
regarding the properties of temporal networks, questions 
which have been addressed through the formulation of 
theoretical models, aimed at explaining both the tempo- 
ral patterns observed and their effects on the correspond- 
ing integrated networks [T7H^ . 



Especially interesting in this perspective is the activity 
driven social network model recently introduced by Perra 
et al. |21| . aimed in particular to capture the relation 
between the dynamical properties of social temporal net- 
works and the topological properties of their correspond- 
ing time projections. The key element in the definition 
of this model is the observation that the formation of so- 
cial interactions is driven by the activity of individuals, 
urging them to interact with their peers, and by the em- 
pirical fact that different individuals show different levels 
of social activity 1211 . Based in the concept of activity po- 
tential, defined as the probability per unit time that an 
individual engages in a social activity, Ref. [H] proposed 
an activity driven social network model, in which indi- 
viduals start interactions, that span for a fixed length of 
time At, with probability proportional to their activity 
potential. The model output is thus given by a sequence 
of graphs, depending on the distribution F{a) of the ac- 
tivity potential, which are updated every time interval 
At. The topological properties of the integrated activity 
driven network were related in Ref. [21 at the level of 
the degree distribution, which, by means of approximate 
arguments, was shown to be proportional to the activity 
potential distribution F{a). However, despite the interest 
of the model, expressions for the rest of topological ob- 
servables have been still lacking, a fact that hampers its 
possible validation as a generator of realistic integrated 
social networks, as well as the identification of the par- 
ticular role that integration time has on the behavior of 
dynamical processes running on top of the temporal net- 
work [221 [Ml. 



Here we address the study of the aggregated network 
generated by the activity driven model, obtained by in- 
tegrating the temporal network up to a given time T, 
by considering a mapping of the integrated network to 
a hidden variables model [21] depending on the activity 
potential distribution and the considered time T. We 
obtain a set of expressions for the degree distribution, 
degree-degree correlations and clustering coefficient of 
the aggregated network that are exact in the limit of 



large network size TV — >■ oo and finite time T, and which 
are amenable to analytic asymptotic expansions in this 
same limit. The expressions obtained, confirmed by nu- 
merical simulations, corroborate the basic assumption of 
the activity driven model linking social activity with net- 
work topology. Moreover, the formalism proposed can be 
extended to generalizations of the activity driven model, 
opening thus the path to the analytical solution of model 
extensions, aimed at better reproducing the topological 
features of real integrated social networks. 

We have organized the paper: Section |III| defines the 
activity driven network model. Sec. |III| gives a brief re- 



view of the hidden variables network formalism; Sec. IV 
deduces the mapping of the activity driven model onto 
this formalism. In Sec. |V] we compute the topological 
properties of the integrated activity network model as a 
function of time T, providing exact expressions as well as 
asymptotic results in the limit T /N small, for different 
forms of the activity potential distribution. In Sec. |Vl] 
we discuss the extension of the formalism presented to 
possible variations of the activity driven model. Finally, 
we summarize our results and conclusions in Sec. IVIIl 



II. THE ACTIVITY DRIVEN NETWORK 
MODEL 

The activity driven network model proposed in 
Ref. |21j is defined in terms of N individuals i (agents), 
each one of them characterized by her activity potential 
fli, defined as the probability that she engages in a so- 
cial act/connection with other agents per unit time. The 
activity of the agents is a (quenched) random variable, 
extracted from the activity potential distribution F{a), 
which can take a priori any form. The model is defined by 
means of a synchronous update scheme, time being mea- 
sured in units of the life span of each connection At. It 
proceeds by creating a succession of instantaneous net- 
works Qt, t — 0, At, 2 At, . . . , nAt, ... At a given time 
t, all previous edges are deleted and we start with N 
disconnected individuals. Each one of them is checked 
and becomes active with probability a^At. Active agents 
generate m links (start m social interactions) that are 
connected to m other agents selected uniformly at ran- 
dom. Finally, time is updated as i — >■ t -|- At. This 
procedure implies that all edges in the temporal network 
have the same constant time duration At. In order to 
avoid complications due to the differences in the number 
of emitted and received connections arising from using a 
synchronous approacnH here we consider a probabilistic 
recipe for the instantaneous network construction: Each 
microscopic time step Ai, we choose N agents, uniformly 
at random, and check sequentially each one of them for 



activation and eventual link emission. We avoid self and 
multiple connections. 

To simplify the analytical calculations performed be- 
low, in the following we choose At = m = 1. Both 
quantities can be however restored by a simple rescaling 
of the activity potential and the integration time T. We 
notice that imposing At = 1 implies restricting the ac- 
tivity potential to be probability, and thus to be limited 
in the interval a e [0, 1]. 



III. HIDDEN VARIABLES FORMALISM: A 
SHORT REVIEW 

The class of network models with hidden variables was 
introduced in Ref. [24, (see also [25l[26]) as a generaliza- 
tion of the random network Gilbert model [2| , in which 
the probability of connecting two vertices is not constant, 
but depends on some intrinsic properties of the respec- 
tive vertices, their so-called hidden variables. This class 
of models is defined as follows: Starting from a set of N 
disconnected vertices and a general hidden variable h, we 
construct an undirected network with no self-edges nor 
multiple connections, by applying these two rules: 

1. To each vertex i, a variable hi is assigned, drawn 
at random from a probability distribution p{h). 

2. For each pair of vertices i and j, i ^ j, with hid- 
den variables hi and hj, respectively, an edge is 
created with probability r{hi,hj), the connection 
probability, which is a symmetric function bounded 
by < r{h,h') < 1. 

Each model in the class is fully defined by the func- 
tions p{h) and r{h,h'), and all its topological proper- 
ties can be derived as a function of these two param- 
eters. These topological properties are encoded in the 
propagator g(k\h), defined as the conditional probabil- 
ity that a vertex with hidden variable h ends up con- 
nected to k other vertices. The propagator is a normal- 
ized function, '^f^g{k\h) = 1, whose generating function 
g{z\h) — ^^ z^g{k\h) fulfills the general equation [23] 

\ng{z\h) = NY,p{h')\n[\ - (1 - z)r{h,h')] . (1) 

h' 

From this propagator, expressions for the topological 
properties of the model can be readily obtained [2H : 



Degree distribution: 



P{k)^Y.9ik\h)p{h). 



(2) 



• Degree correlations, as measured by the average de- 
gree of the neighbors of the vertices of degree k, 
k'^'^ik) [57]: 



^ Indeed, in a synchronous scheme, every time step an agent fires at 
most one connection, but can receive a number n of connections, 
given trivially by a binomial distribution. 



F"(A:) = 1 + 



1 



P(fc) 



Y,pih)g{k\h)k--{h), (3) 



where we have defined 

N 






and 



k{h) ^N^p(h')r{h,h'), 



(5) 



which is the average degree of the vertices with hid- 
den variable h. 

• Average clustering coefficient (c), defined as the 
probability that two vertices are connected, pro- 
vided that they share a common neighbor [55] 



{c) = Y,p{h)c{h), 



(6) 



where we have defined 

-c{h)^ Y,p{h'\h)r{h',h")p{h"\h), (7) 



h'h" 



and 



p{h'\h) 



Np{h')r{h, h') 
Hh) 



(8) 



Additionally, one can define the clustering spec- 
trum, as measured by the average clustering co- 
efficient of the vertices of degree fc, c{k) (571 HH] 



c{k) = j^^Y.P^^)9{k\h)-c{h), 



(9) 



IV. MAPPING THE INTEGRATED NETWORK 
TO A HIDDEN VARIABLES MODEL 

The activity driven network model generates a time se- 
ries of instantaneous sparse networks, with an average de- 
gree {k)t — 2{a), where (a) = J2a^^i^)- The integrated 
network at time T is constructed by performing the union 
of the instantaneous networks, i.e. Qt = ^'tLoGt- In 
this integrated network, vertices i and j will be joined 
by an edge if there has ever been a connection created 
between them in any of the instantaneous networks at 
< t < T. The key point to map the integrated net- 
work to a hidden variables model resides in computing 
the probability nr(i, j) that two vertices i and j become 
eventually joined at time T. This probability is given by 
nT('*, j) = 1 — Qr^hj), where Qrihj), the probability 
that no connection has ever been created between agents 
i and j up to time T, can be calculated as follows: At 
time T, an agent i will have become active z times with 
probability Pt{z). Given the definition of the model, at 
time T we have selected TN agents to check for activa- 
tion. The number of times z that agent i has become 
active will be given by the binomial distribution 



Pt{z) = 



TN 



1- 



N 



TN-z 



(10) 



and analogously for agent j. Now, vertices i and j will 
be connected in the integrated network if at least one of 
the links generated from i reaches j, or vice- versa. Since 
every time that she becomes active, an agent creates a 
connection targeted to a randomly chosen peer, the prob- 
ability Qrihj) is given by 



MhJ)^J2PT{z^)PT{z,)(l 






1 

N 

TN 



1 



1 

N 



(11) 



where we have performed the summation using the prob- 
ability distribution in Eq. (10). We see now that the 
probability that agents i and j are connected in the in- 
tegrated network at time T depends only on their re- 
spective activity potentials Oi and a^, which are random 
variables with distribution F{a). The mapping to a hid- 
den variables network is thus transparent: 



Hidden variable: h — > 



a. 



• Distribution of hidden variables: p{h) — > F{a). 

• Connection probability: r{h,h') — > nT(a, a'). 

At very large times, the integrated network emerging 
from the activity driven model will trivially tend to a 
fully connected network. Interesting topology will thus 
be restricted to the limit of small T compared with the 



network size N . In this limit, Eq. (11) can be simplified, 
yielding 



nT(a,a') = 1 - QT{a,a) ~ 1 - 
~ 1 — exp [— A(a + a')] , 



1 



(a + a') 

iV2 



TN 



(12) 



where we have neglected terms of order 0{N '^) and de- 
fined the parameter 



A = 



N' 



(13) 



An explicit calculation of the connection probability for 
a factor m > 1 and a time interval At ^ 1 can be easily 
performed; in the limit of large N and constant A, the 
only change ensuing is a rescaling of time, T — > Ttti, the 
value of At becoming canceled in the process of taking 
the hmit A — ^ 0. 



V. TOPOLOGICAL PROPERTIES OF THE 
INTEGRATED ACTIVITY DRIVEN NETWORK 

Here we will apply the formalism presented in Sec. |III| 
to provide analytic expressions characterizing the topol- 
ogy of the integrated network resulting from the activity 
driven model. For the sake of concreteness, we will focus 
in the following activity potential distributions, in the 
continuous a limit: 



• Constant activity: 

F{a) = Sa.acn ^^^^ < Oq < 1. 

• Homogeneous activity: 

F{a) = l/flmax, with < a < amax < 1- (14) 

• Power-law distributed activity: 

F{a)^{"/-l)e^-^a-'', withae[£,l]. (15) 

In the last case, where we consider 7 > 2, in accordance 
with experimental evidence |21) . we have introduced a 
lower cut-off < e ^ 1 in order to avoid dangerous 
divergences in the vicinity of zero. 

A. Degree distribution 

In order to compute the degree distribution, we have to 
solve and invert the generating function equation Eq. (nj), 
an almost impossible task to perform exactly, except in 
the case of very simple forms of the activity potential 
distribution. So, in the case of constant activity, F{a) — 
Sa,aoi we have 



g{z\ao) = [zTlriaQ, oo) + (1 - T^riao, ao)) 



N 



(16) 



which corresponds to the generating function of a bino- 
mial distribution ,30j . Therefore, in the limit of large N 
and constant A, the degree distribution takes the Poisson 
form 



Prik) 



fc! 



(17) 



with parameter fi — N (l ~ e^^'*"^°), which, for fixed T 
and large N, can be approximated as /i ~ 2TaQ. 

For a nontrivial activity distribution F(a), we must re- 
sort to approximations. We therefore focus in the inter- 
esting limit of small A, which corresponds to fixed T and 
large N, which is the one yielding a non-trivial topol- 
ogy. In this limit, we can approximate the connection 
probability as 



11^(0, a ) ~ A(a + a ). 



(18) 



Introducing this expression into Eq. ([I]) and performing 
a new expansion at first order in A, we obtain [24j 

\ng{z\a) ~ (1 - z)AiV^ F(a')(a + a') (19) 



= {l-z)XN{a+{a)). 



(20) 



The generating function of the propagator is a pure ex- 
ponential, which indicates that the propagator itself is a 
Poisson distribution RHI, i-e. 



ma) 



-T(a+(a)) F(« +(«))]' 



r(A: + l) 



(21) 



where r(a;) is the Gamma (factorial) function [21]. From 
Eq. (pi) we obtain the general expression for the degree 
distribution 



Prik) 



rpk 



nk+i)^ 



^f(a)[a+(a)]"e-^(''+<'^». 



(22) 



In the case of a homogeneous activity distribution, 
F(a) = ttmax"^, for which (a) — a,„ax/2, we can inte- 
grate directly Eq. ( 22 1 , to obtain 



p .,^ r(fc + i,r(a))-r(fc + i,3r(a)) 

where r(x, z) is the incomplete Gamma function |31j . 

More complex forms of the activity distribution do not 
easily yield to an exact integration, and more approxima- 
tions must be performed. In particular, the asymptotic 
form of the degree distribution can be obtained by per- 
forming a steepest descent approximation. Thus, we can 
write 



PT{k) = 



1 



r(fc + i) 

where we have defined 



F(a)e'^(°)da, 



(j){a) = k\n[T{a+{a))]-T{a + {a)). 



(24) 



(25) 



The function 0(a) has a sharp maximum around aM = 
^ — (a). Performing a Taylor expansion up to second 

order, we can write (j){a) ~ 4>{aM) ~ ^\a ~ om ] , with 
4>{aM) == fcln(fc) - k. Now, for T^/fc > 1, the tunc- 

T^ r ]2 

tion e^ ^T^i^-^Mi is strongly peaked around the maxi- 
mum aM] therefore we can substitute the activity poten- 
tial by its value at the maximum, to obtain 



PT{k) 



e't>('"'')F{aM) 
r(fc + l) 

^/2^,k^e-^ 
Tr(fc + 1) 



^-^[a-aM\ ^^ 



00 

k 



(26) 



where we have extended the integration limits to plus and 
minus infinity. In the large k limit, we can use Stirling's 



approximation, T(k 
asymptotic form 



Prik) - ^F [^ 



to obtain the 



(27) 



In this expression we recover, using more rigorous ar- 
guments, the asymptotic form of the integrated degree 
distribution obtained in Ref. [21]. The limits of valid- 
ity of this expression are however now transparent, being 
exphcitly iV > T > 1 and T^ > fc > 1. 

For the case of constant activity, F(a) = Sa,ao: 
the asymptotic form of the degree distribution is 
Prik) ^ Sk^Tao/T, while the exact form is a Poisson dis- 
tribution centered at 2Tao. For a uniform activity, on 
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FIG. 1. (color online) Rescaled degree distribution Prik) for 
integrated networks corresponding to different values of T, 
with power law activity distribution with exponents 7 = 3 
and 2.25. Network size N = 10^. The behavior predicted by 



Eq. (271 is represented as dashed lines. 



the other hand, the asymptotic prediction is a flat dis- 
tribution, while the exact expression can be quite dif- 
ferent, in particular for large and small values of fc, see 
Eq. (23). For the case of a power-law distributed activ- 
ity, in Fig. [1] we plot the degree distribution Prik) of 
the aggregated network at different values of T for net- 
works of size iV = 10^ and two different values of 7. As 
we can see, for such large networks sizes and values of 
A ~ 10^^ — 10^"^, the asymptotic expression Eq. (27) 



represents a very good approximation to the model be- 
havior. In Fig. [2] we plot the degree distribution for a 
smaller network size N = 10'^. As one can see, a numer- 



ical integration of Eq. ( 22 ) recovers exactly the behavior 



of Prik) even for small values of k. With such small net- 
work size, however, the asymptotic prediction of Eq. (27) 
is less good, as shown in the inset of Fig. [2] 



FIG. 2. (color online) Degree distribution Prik) for inte- 
grated networks corresponding to different values of T, with 
power law activity distribution with exponent 7 = 2.25. Net- 
work size A'^ = 10^. The result of a numerical integration of 



Eq. (22 1 is showed as c ontinuous lines. Inset: Rescaled Prik) 
shown against Eq. (271, dashed in blue. 



a constant activity potential, F(a) = Sa^aa-, we have 
P-r(fc) = 5(A:|ao). Since in this case *(A) = e"^°° 
have 



we 



fc«"(fc) = l + Ar[i 



-2Aao 



] ~ 1 -f 2rao 



(31) 



where the last expression corresponds to the limit of small 
A. This function is independent of /c, indicating that 
the integrated network corresponding to constant activ- 
ity potential has no degree correlations. 

For more complex forms of F(a) , we resort to an expan- 
sion in powers of A to obtain an approximate expression, 
which at lowest order takes the form 



A:?;"(a) 



AA^ 



{{c?) + {af + 2a(a)] . (32) 



B. Degree correlations 

We start from Eq. (Isl), which takes the form, as a func- 
tion of time 

^r(a) = A^[l-e"^°*(A)], (28) 

where ^(A) is the Laplace transform 

vI/(A)^^na)e-^^ (29) 

a 

We can now use Eq. ([4| , which leads to the exact expres- 
sion 

trW^^{l-*(A) 7_-,^^^Jr }. (30) 

In order to obtain an explicit expression for fc™(fc) we 
must perform the integral in Eq. ([3]). In the case of 



Inserting this expression into Eq. ([3]) , and considering the 
Poisson form of the propagator Eq. (l2l|) , we can write 



^T{k) ^ 1 + ^\pI^^ ^ ^ I daF{a)g{k - l|a) + 



2T2(a) 



kP{k) 

.2^(^-1 



daaF{a)g{k — l\a) 



1 + T^ 



kP{k) 



at + 2(a) 



T 



where in the last expression we have performed the steep- 
est descent approximation used to obtain Eq. (|26]), and 
a\ — (a^) — (a)^ is the variance of the activity potential 
F{a). In the limit of large k, where P{k — 1)/P{k) ~ 1, 
we have the general form for the degree correlations 



hnn 



(fc)-l 



2{a)+al 



T 



(33) 
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FIG. 3. (color online) Main: Rescaled average degree of the 
nearest neighbors of the vertices of degree k, kj-^ik), for the 
integrated network with size A*' — 10'' and power law activity 
distribution with 7 — 2.5, for different values of T. The pre- 
diction of Eq. ( 33 1 is dashed in blue. Inset: Average degree of 
the neighbors of the vertices with activity a, k^"{a), for the 



same integrated network. The predictions from Eq. ( 32 1 are 
shown as continuous lines. 



This expression recovers in a natural way the exact re- 
sult for constant activity potential, where a^ = 0. From 
Eq. (33) we conclude that, in general, for an non-constant 



activity distribution, the integrated networks resulting 
from the activity driven model show disassortative mix- 
ing by degree [35], with a fc™(fc) function decreasing as a 
function of k. This disassortative behavior, which can be 
however quite mild in the case of small variance a a , as in 
the case of a power law distributed activity with small e, 
is in any case at odds with the assortative form observed 
for degree correlations in real social networks [1]. 

In Fig. [3] we check the validity of Eq. ( 32 1 and the 



asymptotic form Eq. ( 33 ) in the case of power law dis- 
tributed activity. We observe that the prediction of 
Eq. (32) recovers exactly the model behavior, also in 
the case of small activity a (shown in the inset). The 
degree correlation, kip{k)^ also correctly captured by 



the asymptotic form Eq. ( |33[ ). Note however that, since 
the variance da is small (of order £'^~^ for 7 < 3 and or- 
der e^ for 7 > 3), the net change in the average degree of 
the neighbors is quite small, and the integrated network 
can be considered as approximately uncorrelated without 
incurring in a gross error. 



C. Clustering coefficient 

The expression of the clustering spectrum at time T, 
cr(fc), takes the form, from Eq. (jof 

CT(fc) = -^ Y. F{o)9(k\a)cT{<^). (34) 



FIG. 4. (color online) Main: Rescaled clustering coefficient 
of the nodes of degree k, CT{k), of the integrated network 
with size A^ = 10'' and power law activity distribution with 
7 — 2.5, for different values of T. The prediction of Eq. (391 
is dashed in blue. Inset: Clustering coefficient of the nodes 
with activity a, cria), of the same integrated network. The 



predictions from Eq. ( 37 1 are shown as continuous lines. 



Using Eqs. (I?]), (IS]) and the expression for fc(a), we can 
write the exact form 



cria) = 1 - 



*(A) - e-^''^(2A) 
l-e-^'^*(A) 



(35) 



Again in the simplest case of a constant activity poten- 
tial, F{a) — (5a,aoi we have cria) = 1 — e~^^°-°, which 
leads to a clustering spectrum 

2Tao 



CT{k) = {c)t = 1- e 



-2Aao 



N 



(36) 



where the last expression is valid for small A. The clus- 
tering spectrum is in this case constant, and equal to 
the average clustering coefficient. For fixed time T, it is 
inversely proportional to the network size, in correspon- 
dence to a purely random network. It increases with T, 
saturating at (c)oo ~ 1 for a fully connected network in 
the infinite time limit. 

For a general activity potential distribution, we need 
to perform again an expansion in A, which in this case 
takes the form, at first order in A, 

o\ 

c{a)^ ^[{a^) + a{a)]. (37) 



Inserting this form into Eq. (34), and performing the 



same steepest descend approximation applied in Eq. ( 33 1 , 
we obtain 



m 



2r2 p{k - 1) 

~n' kP{k) 



{{a' 



') + («) 



(38) 



In the limit of large fc, we obtain the general form of the 
clustering spectrum, valid for any activity potential. 



T 



N 



N 



(39) 



In Fig. [4] we plot the clustering coefficient as a function 
of the degree (main) and the activity (inset), in the case 
of power law distributed activity. We observe that both 
Eq. (37) and Eq. (39) recover correctly the clustering 



coefficient behavior. 



appropriate expression for the attractiveness distribution 

Fib). 



VII. CONCLUSIONS 



VI. MODEL EXTENSIONS 

The potency of the hidden variables formalism we have 
introduced above to solve the activity driven model al- 
lows to easily extended it to tackle the analysis of gen- 
eralized models inspired in the same principles. We can 
consider, indeed, different rules for activation and recep- 
tion of connections. The only limitation to be imposed in 
order to properly implement the formalism is that con- 
nection rules must be local, i.e. involving only proper- 
ties of the emitting and receiving agents. As a simple 
example, we consider a sort of "inverse" activity driven 
model, in which every agent i becomes active with the 
same constant probability a^ = oq and, when active, she 
sends a connection to another agent j, chosen at ran- 
dom with probability proportional to some (quenched) 
random quantity bj, the attractiveness of the node, i.e. 
with probability bj/{b)N. In this case, one can easily re- 
peat the steps of the mapping presented in Sec. |IV| The 
number of times z that agent i becomes active is now 



P^{z) 



TN\ 
z J KN 



/ao 



N 



TN- 



(40) 



and the probability that i and j never become connected 
up to time T is 



Q'T{hJ) = J2PTiz^)PTizJ)l^ 



{b)N 



{b)N 



{b)Ny \ {b)N^J 
exp[-X'{bi + bj)] , 



TN 



The activity driven model represents an interesting ap- 
proximation to temporal networks, providing an prelim- 
inary explanation of the origin of the degree distribution 
of integrated social networks, in terms of the heterogene- 
ity of the agents' activity, and the distribution of this 
quantity. Here we have explored the full relation between 
topology and activity distribution, obtaining analytical 
expressions for several topological properties of the in- 
tegrated social networks for a general activity potential, 
in the thermodynamic limit of large number of agents, 
N —)' 00, and finite integration time T. To tackle this 
issue, we have applied the hidden variables formalism, 
by mapping the aggregated network to a model in which 
the probability of connecting two nodes depends on the 
hidden variable (in this case represented by the activity 
potential) of those nodes. Our analysis is complemented 
by numerical simulations in order to check theoretical 
predictions against concrete examples of activity poten- 
tial distributions. Using our formalism, we can demon- 
strate rigorously that the integrated degree distribution 
at time T takes the same functional form as the activ- 
ity potential distribution, as a function of the rescaled 
degree k/T — (a). This is however an asymptotic result, 
which is well fulfilled for an activity potential power-law 
distributed, as empirically measured in a wide range of 
social interaction settings, which fails for simple constant 
or homogeneous distributions. We also show that the ag- 
gregated networks show in general disassortative degree 
correlations, at odds with the assortative mixing revealed 
in real social networks. The clustering coefficient is low, 
(c) ^ T/N, comparable with a random network. 

Our study opens interesting direction for future work, 
concerning, for example, the clarification of role of inte- 
gration time in the properties of dynamical process on 
activity driven networks, and the possible modifications 
of the activity driven network model, in order to incor- 
porate some properties of real social networks currently 
missed, such as a high clustering coefficient, assortative 
mixing by degree or a community structure [1]. 



where we have defined the new parameter A' = a^T / {b)N 
and, in the last step of the previous expressions, we have 
performed and expansion for large TV a finite T. From 



here. 



we obtain n^(i,j) = 1 — Q^(?,j). As we can see. 



this modified model is exactly mappable to the activity 
driven model (see Eq. (12), with the simple translation 
A -> A'; all the general expression derived above hold thus 
in this case, and can be worked out, upon providing the 
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